5*4 Symmetric Microscopic Model for Iron-Based Superconductors 
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Although iron-based superconductors are multi-orbital systems with complicated band structures, 
we demonstrate that the low energy physics which is responsible for high-Tc superconductivity is 
essentially governed by one effective Hamiltonianwith two almost decoupled orbitals near half filling. 
This underlining electronic structure is protected by the Si symmetry. With repulsive or strong next 
nearest neighbor antiferromagnetic exchange interactions, each single-orbital effective Hamiltonian 
results in a robust Aig s-wave pairing which can be exactly mapped to the d-wave pairing observed 
in cuprates. The classification of the superconducting(SC) states according to the 54 symmetry 
leads to a natural prediction of the existence of two different phases named A and B phases. In the 
B phase, the superconducting order has an overall sign change along c-axis between the top and 
bottom As(Se) planes in a single Fe-(As)Se trilayer structure, which is an analogy of the sign change 
under the 90° degree rotation in the d-wave SC state of cuprates. Our derivation provides a unified 
understanding of iron-pnictides and iron-chalcogenides, and suggests that cuprates and iron-based 
superconductors share identical high-Tc superconducting mechanism. 



I. INTRODUCTION 

Since the discovery of iron-based superconductorsir— , 
there has been considerable controversy over the choice 
of the appropriate microscopic Hamiltonian^i^. The ma- 
jor reason behind such a controversy is the complicated 
multi d-orbital electronic structure of the materials. Al- 
though the electronic structure has been modeled by us- 
ing different numbers of orbitals, ranging from minimum 
two orbitals-, three orbitalsS, to all five d orbital s^' a 
general perception has been that any microscopic model 
composed of less than all five c?-orbitals and ten bands 
is insufficient^. Such a perception has blocked the path 
to understand the superconducting mechanism because 
of the difficulty in identifying the key physics responsible 
for high Tc- Realistically, in a model with five orbitals, 
it is very difficult for any theoretical calculation to make 
meaningful predictions in a controllable manner. 

Iron-based superconductors include two families, iron- 
pnictidesi^ and iron-chalcogenides^. They share many 
intriguing common properties. They both have the high- 
est TcS around SOK^i^iii"— . The superconducting gaps 
are close to isotropic around Fermi surfacesi^"— and the 
ratio between the gap and Tc, 2A/Tc, are much larger 
than the BCS ratio, 3.52, in both families. However, the 
electronic structures in the two families, in particular, 
the Fermi surface topologies, are quite different in the 
materials reaching high Tc. The hole pockets are absent 
in iron-chalcogenides but present in iron-pnictidesi^^— . 
The presence of the hole pockets has been a necessity for 
superconductivity in the majority of studies and mod- 
els which deeply depend on the properties of Fermi sur- 
faces. Therefore, the absence of the hole pockets in iron- 
chalcogenides causes a strong debate over whether both 
families belong to the same category that shares a com- 
mon superconducting mechanism. Without a clear mi- 
croscopic picture of the underlining electronic structure, 
such a debate can not be settled. 



Observed by angle-resolved photoemission microscopy 
(ARPES), a very intriguing property in the SC states 
of iron-pnictides is that the SC gaps on different Fermi 
surfaces are nearly proportional to a simple form factor 
coskxCosky in reciprocal space. This form factor has been 
observed in both 123iL22,-^ and 111^^1^ families of iron- 
pnictides. Just like the d-wave form factor cosk^ — cosky 
in cuprates, such a form factor indicates that the pair- 
ing between two next nearest neighbour iron sites in real 
space dominates. In a multi orbital model, many theo- 
retical calculations based on weak coupling approaches 
have shown that the gap functions are very sensitive to 
detailed band structures and vary significantly when the 
doping change o^'^^" — . The robustness of the form fac- 
tor has been argued to favor strong coupling approaches 
which emphasize electron-electron correlation or the ef- 
fective next nearest neighbour (NNN) antiferromagnetic 
(AF) exchange coupling J^^^— as a primary source of 
the pairing force. However, realistically, it is very dif- 
ficult to imagine such a local exchange interaction re- 
mains identical between all c?-orbital electrons if a multi 
d-orbital model is considered. 

In this paper, we demonstrate that the underlining 
electronic structure in iron-based superconductors, which 
is responsible for superconductivity at low energy, is es- 
sentially governed by a two orbital model obeying the S4 
symmetry. The two orbital model includes two nearly de- 
coupled single-orbital parts that can be mapped to each 
other under the S4 transformation. This electronic struc- 
ture stems from the fact that the dynamics of dxz and 
dyz orbitals are divided into two weakly coupled groups 
that are separately coupled to the top and bottom As(Se) 
planes in a in a single Fe-(As)Se trilayer structure. The 
two groups can thus be treated as a iso-spin. 

The underlining electronic structure becomes transpar- 
ent after performing a gauge mapping in the five orbital 
modeli^. The gauge mapping also reveals the equivalence 
between the Aig s-wave pairing and the d-wave pairing. 
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After the gauge mapping, the band structure for each 
^4 iso-spin component is characterized by Fermi surfaces 
located around the anti d-wave nodal points in Brillouin 
zone, corresponding to the sublattice periodicity of the 
bipartite iron square lattice as shown in Fig [if a). In the 
presence of an AF exchange coupling J2 or an effective 
on-site Hubbard interaction, the d-wave pairing defined 
in the sublattices can be argued to be favored, just like 
the case in cuprates. The d-wave pairing symmetry maps 
reversely to a Aig s-wave pairing in the original gauge 
setting. These results provide a unified microscopic un- 
derstanding of iron-pnictides and iron-chalcogenides and 
explain why an s-wave SC state without the sign change 
on Fermi surfaces in iron chalcogenides driven by repul- 
sive interaction can be so robust. More intriguingly, since 
the different gauge settings do not alter any physical mea- 
surements, the results suggest that in the Aig s-wave 
state, for each 6*4 iso-spin component, there is a hid- 
den sign change between the top As(Se) and the bottom 
As(Se) planes along c-axis. 

The 5*4 symmetry adds a new symmetry classifica- 
tion to the SC states. For example, even in the Aig 
s-wave pairing state, there are two different phases called 
A and B phases, with respect to the 5*4 symmetry. In 
the A phase, the relative SC phase between the two 54 
iso-spin components is zero while in the B phase, it is 
TT. Therefore, there is an overall tt phase shift between 
the top As(Se) and the bottom As(Se) planes in the B 
phase along the c-axis. Such a sign change should be de- 
tectable experimentally. This property makes iron-based 
superconductors useful in many SC device applications. 
An experimental setup, similar to those for determining 
the d-wave pairing in cuprates^"—, is proposed to de- 
tect the TT phase shift. The detection of the sign change 
will strongly support that cuprates and iron-based super- 
conductors share identical microscopic superconducting 
mechanism. 

The paper is organized in the following way. In Sec- 
tion II, we perform a gauge mapping and discuss the 
emergence of the underlining electronic structure. In Sec- 
tion III, we show that the underlining electronic struc- 
ture can be constructed by a two orbital model obeying 
the Si symmetry and discuss many general properties in 
the model. In Section IV, we discuss the classification 
of the SC states under the S4 symmetry and propose a 
measurement to detect the tt phase shift along c-axis be- 
tween the top and bottom As(se) planes. In Section V, 
we discuss the analogy between iron-based superconduc- 
tors and cuprates. 



II. GAUGE MAPPING AND THE 
EQUIVALENCE OF S-WAVE AND d-WAVE 
PAIRING 

Gauge Mapping: We start to ask whether there is an 
unidentified important electronic structure in iron-based 
superconductors in a different gauge setting. Giving a 




FIG. 1. (a) the square lattice structure of a single iron layer; 
one cell includes two Fe ions identified with different filled 
black balls which form two sublattices. We use x — y coor- 
dinate to mark the original tetragonal lattices and x' — y' to 
mark the sublattice direction, (b) the gauge transformation is 
illustrated. The balls with red circles are affected by the gauge 
transformation, (c) and (d) the mapping from the s-wave to 
the d-wave pairing symmetry by the gauge transformation. 

translational invariant Hamiltonian that describes the 
electronic band structure of a Fe square lattice, 

-^0 = ^ U],al3fta,af]P,<y, (1) 

where z,j label Fe sites, a,/3 label orbitals and a labels 
spin. We consider the following gauge transformation. 
As shown in Figma,b) , we group four neighbouring iron 
sites to form a super site and mark half super sites by red 
color. The gauge transformation, [/, adds a minus sign 
to all Fermionic operators fia.a at every site i marked 
by red color. After the transformation, the Hamiltonian 
becomes 

id^ = U+HqU. (2) 

The gauge mapping operator U is an unitary operator 
so that the eigenvalues of Tdo are not changed after the 
gauge transformation. It is also important to notice 
that the mapping does not change standard interaction 
terms, such as conventional electron-electron interactions 
and spin-spin exchange couplings. Namely, for a general 
Hamiltonian including interaction terms _/d/, under the 
mapping, 

H = Ha + Hi ^ H' = U+HU = Hq + Hi- (3) 

It is also easy to see that every unit cell of the lat- 
tice in the new gauge setting has four iron sites. The 
original translational invariance of a Fe-As(Se) layer has 
two Fe sites per unit cell. As we will show in the fol- 
lowing section, the doubling of the unit cell matches the 
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FIG. 2. (Three^, fiv©i2) orbital models: (a,e) the Fermi 
surfaces, (b,f) the band dispersion along the high symmetry 
lines, (c,g) the Fermi surfaces after the gauge transformation, 
(d,h) the band dispersions along the high symmetry lines after 
the gauge transformation. The hopping parameters can be 
found in the above two references. 



true hidden unit cell in the electronic structure when the 
orbital degree of freedom is considered. This is the fun- 
damental reason that the new gauge happens to reveal 
the underlining electronic structure. 

Equivalence of s-wave and d-wave pairing: The gauge 
mapping has another important property. As shown 
in Figlljcjd), this transformation maps the Aig s-wave 
cos{kx)cos{ky) pairing symmetry in the original Fe lattice 
to a familiar d-wave cosk'^ — cosk'y pairing symmetry de- 
fined in the two sublattices, where (fc^, ky) and (fc^, k'y) la- 
bel momentum in Brillouin zones of the origin lattice and 
sublattice respectively. A similar mapping has been dis- 
cussed in the study of a two-orbital iron ladder model^i^ 
to address the equivalence of s-wave and d-wave pairing 
symmetry in one dimension. 

In an earlier paper-2^, one of us and his collabora- 
tor suggested a phenomenological necessity for achieving 
high Tc and selecting pairing symmetries: when the pair- 
ing is driven by a local AF exchange coupling, the pairing 
form factor has to match the Fermi surface topology in 
reciprocal space. If this rule is valid and the iron-based 
superconductors are in the Aig s-wave state, we expect 
that the Fermi surfaces after the gauge mapping should 
be located in the d-wave anti-nodal points in the sublat- 
tice Brillouin zone. 

Band structures after gauge mapping: There have been 
various tight binding models to represent the band struc- 
ture of Hq. In FiglU we plot the band structure of iJo 
and the corresponding //q for two different models: a 
maximum five-orbital model for iron-pnictides^", and a 
three-orbital model constructed for electron-overdoped 
iron-chalcogenides^^. 

As shown in Fig (21 although there are subtle differ- 
ences among the band structures of H'q , striking common 
features are revealed for both models. First, exactly as 
expected, all Fermi surfaces after the gauge mapping are 
relocated around X', the anti-nodal points in a standard 
d-wave superconducting state in the sublattice Brillouin 



zone. This is remarkable because a robust d-wave super- 
conducting state can be argued to be favored in such a 
Fermi surface topology in the presence of repulsive in- 
teraction or nearest neighbour (NN) AF coupling in the 
sublattice^i^. If we reversely map to the original gauge, 
the original Hamiltonian must have a robust s-wave pair- 
ing symmetry. Therefore, an equivalence between the 
Aig s-wave and the d-wave pairing is clearly established 
by the gauge mapping. 

Second, the bands previously located at the different 
places on the Fermi surface are magically linked in the 
new gauge setting. In particular, the two bands that con- 
tribute to electron pockets are nearly degenerate and in 
the five orbital model, the bands that contribute to hole 
pockets are remarkably connected to them. Considering 
the fact that the unit cell has four iron sites in the new 
gauge setting, this unexpected connections lead us to be- 
lieve that in the original gauge, there should be just two 
orbitals which form bands that make connections from 
lower energy bands to higher energy ones and determine 
Fermi surfaces. Moreover, the two orbitals should form 
two groups which provide two nearly degenerate band 
structures. Finally, since the mapping does not change 
electron density, Figl^] reveals the doping level in each 
structure should be close to half filling. 

In summary, the gauge mapping reveals that the low 
energy physics is controlled by a two orbital model that 
produces two nearly degenerated bands. 



III. THE CONSTRUCTION OF A 
TWO-ORBITAL MODEL WITH THE 54 
SYMMETRY 

With above observations, we move to construct an ef- 
fective two orbital model to capture the underlining elec- 
tronic structure revealed by the gauge mapping. 

Physical picture: Our construction is guided by the fol- 
lowing several facts. First, the d-orbitals that form the 
bands near the Fermi surfaces are strongly hybridized 
with the p-orbitals of As(Se). Since the d^'z and dyiz 
have the largest overlap with the px' and Pyi orbitals, it 
is natural for us to use d^'z and dy'z to construct the 
model. Second, in the previous construction of a two- 
orbital model, the C/^y symmetry was used^. The C^y 
symmetry is not a correct symmetry if the hopping pa- 
rameters are generated through the p-orbitals of As{Se). 
Considering the As(Se) environment, a correct symme- 
try for the d-orbitals at the iron-sites is the S'4 symmetry 
group. Third, there are two As(Se) planes which are sep- 
arated in space along c-axis. Since there is little coupling 
between the p orbitals of the two planes and the hop- 
pings through the p-orbitals are expected to dominate 
over the direct exchange hoppings between the d-orbitals 
themselves, the two orbital model essentially could be 
decoupled into two weakly coupled one orbital models. 
Finally, the model should have a translational invariance 
with respect to the As(Se) plane. 
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FIG. 3. A sketch of the d^/^ and dyi^ orbitals, their orien- 
tations and their couphng into the two As(Se) layers. The 
hopping parameters are indicated: the nearest neighbor hop- 
ping is marked by tix,iy, the next nearest neighbor hoppings 
are t2 and t'2 due to the broken symmetry along two differ- 
ent diagonal directions, the third NN hopping is marked by 
t3x,3y The coupling between two layers is marked by the 
nearest neighbor hopping tc- 



With above guidelines, it is very natural for us to di- 
vide the two d-orbitals into two groups as shown in Fig|31 
One group includes the d^'z in the A sublattice and the 
dy'z in the B sublattice, and the other includes the d^'z in 
the B sublattice and the dy'z in the A sublattice, where A 



and B label the two sublattices of the iron square lattice 
as shown in FiglIJa). The first group strongly couples 
to the p-orbitals in the up As(Se) layer and the second 
group couples to those in the bottom As(Se) layer. We 
denote Cia and dia as Fermionic operators for the two 
groups respectively at each iron site. 

5*4 symmetry and the two-orbital model: Without turn- 
ing on couplings between the two groups, we seek a gen- 
eral tight binding model to describe the band structure 
based on the Si symmetry. The S'4 transformation maps 
Cia- to dicr. If we define the corresponding operators in 
momentum space as Cka and dka, the S'4 transformation 
takes 



dka 



Ck'+Qa 



(4) 



where k' ~ {ky,—kx) and Q — {tt,7t) for given k — 
{kx, ky ). 

Now, we consider a tight binding model for the first 
group. Here we limit the hopping parameters up to 
the third NN (TNN). As illustrated in FigH the tight 
binding model can be approximated by including NN 
hoppings, tlx, tiy, NNN hoppings, t2, t'2, and TNN 
hoppings, t^x and t^y. The longer range hoppings can 
be included if needed. For convenience, we can define 
tis — {tlx + tiy)/2, tid ~ {tlx — *ii/)/2, t2s = {t2 + ^2)72 
and t2d = {t2 - t'2)/'^, hs = {tsx + hy)/'^ and t^d = 
{tzx — ^3i/)/2, where the labels, s and d, indicate s-wave 
(hopping parameter is symmetric under the 90° degree 
rotation) and d-wave (hopping parameter changes sign 
under the 90° degree rotation) type hoppings respec- 
tively. A general tight binding model can be written 
as 



J 



Ho^one = ^ 2[tis{cOskx + COsky) - ^ + tid{cOskx - COsky)]c^^Cka + 4[t2sCOskxCOskyC^^Cka + t2dSinkxSinkyC^^Ck+Qa] 



'2[t'is{cos2kx + cos2ky) + t3d{cos2kx - cos2ky)]c^^Cko 



(5) 



We can apply the S'4 transformation on Hq to obtain the 
tight binding model for the second group. It is straight- 
forward to see that the transformation invariance requires 
tis,t2d and t^d to change signs. Therefore, the two-orbital 
model is described by 



Ho^two = 2_^['it2sCOskxCOsky - fi]{c+^Ckcr + dl^dka) 
k(T 

+2tis{coskx + cosky){cl^Ckc, - d^^dka) 



+2tid{coskx ~ cosky){cl^Cka + dj^dka) 

+At2dSinkxSinkyic'^^Ck+Qa- - dta'^k+Qa) 

+2t3s{cos2kx + C0s2ky){cl^Cka + d'l^dka) 

+2t3d{cos2kx - COs2ky){cl^Ckcr - d'^^dka) 
+... (6) 

Ho,two provides the following energy dispersions for the 
two orbitals. 



= ek ± 2t3d{cos2kx — cos2ky) + 4-\/ t^^siri^xsirfiy + [ 



tis{coskx + cosky) ± tid{coskx — cosky) 



Eh± = Cfc ± 2t3d{cos2kx — cos2ky) — 4y t^^siii^xsiri^y + [ 



.tis{coskx + cosky) ± tid{coskx — cosky) 



y']2 



(7) 



(8) 
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FIG. 4. The Fermi surfaces of each component when only 
parameters tis,t2d and t2s are considered. The layout exactly 
follows Fig[2] The parameters are ti = 0.24, t2 = 0.52 and 
jj, — —0.273. The only different parameter between (a) and 
(e) is t'2 with t'2 = -0.1 in (a) and t'2 = -0.2 in (e). 



where = 'it2sCoskxCosky + 2t^s{cos2kx + cos2ky) — fi. 

Ee± captures the electron pockets at M points and 
Ei^± captures the hole pockets at F points. Based on the 
previous physical picture, iis, t2s and t2d should be the 
largest parameters because they are generated through 
the p-orbitals. In FigHl we show that by just keeping 
these three parameters, the model is already good enough 
to capture the main characters of the bands contribut- 
ing to Fermi surfaces in the five-orbital model. After 
performing the same gauge mapping, this Hamiltonian, 
as expected, provides pockets located at X' as shown in 
FigH 

Now we can turn on the couplings between the two 
groups. It is straightforward to show that the leading 
order of the couplings that satisfies the S4 symmetry is 
given by 

Ho,c = ^ 2tc{coskx + cosky){cl^dka + h.c). (9) 

k 

This term can be generated through the NN hopping 
between two orbitals with the same orbital characters. 
Due to the form factor (coskx -\- cosky), the coupling has 
little effect on the electron pockets while has visible effect 
on the hole pockets. Without this coupling, the two hole 
pockets are not only degenerate at F but also have four 
degenerate points along the two T — X directions. The 
coupling can lift the later degeneracy on the hole pockets. 

Combining the Hq^iwo and -ffo.c we obtain an effective 
Si symmetric two orbital model whose band structure is 
described by 



0,eff 



O,two 



O.c. 



(10) 



The c and d Fermionic operators can be viewed as two 
iso-spin components of the Si symmetry. 

General properties of the model: The above model is 
good enough to quantitatively describe the experimen- 
tal results measured by ARPESili^i^^i^. Although the 
hopping parameters are dominated by tis,t2d and i2s, 



other parameters can not be ignored. For example, at 
the same M points, there is energy splitting between two 
components, which indicates the existence of a sizable 
tid- To match the detailed dispersion of the bands, the 
TNN hoppings have to be included. The existence of 
the TNN hoppings may also provide a microscopic jus- 
tification for the presence of the significant TNN AFM 
exchange coupling J3 measured by neutron scattering in 
iron-chalcogenides^^^i^. 

While the detailed quantitative results for different 
families of iron-based superconductors will be present 
elsewhere^, we plot a typical case for iron-pnicitides 
with parameters tis = 0.4, = — 0.03,t2s — 0.3,i2(i = 
0.6, hs = 0.05, tzd = -0.05 and ^ = -0.3 in figl^a-d). In 
Figini;a,b), the coupling U = 0. In Figini;c,d), tc = 0.02. 
It is clear that the degeneracy at the hole pockets along 
T — X direction is lifted by tc- The Fermi surfaces in FiglS] 
are very close to those in the five orbital model"'^^. This 
result really suggests the coupling between two orbitals 
is rather small. 

There are several interesting properties in the model. 
First, it is interesting to point out that we can make 
an exact analogy between the 6*4 transformation on its 
two iso-spin components and the time reversal symme- 
try transformation on a real 1/2-spin because Si = —1. 
This analogy suggests in this S'4 symmetric model, the 
degeneracy at high symmetric points in Brillouin zone is 
the type of the Kramers degeneracy. 

Second, in this model, if the orbital degree of freedom 
is included, the true unit cell for each iso-spin component 
includes four irons. The gauge mapping in the previous 
section exactly takes a unit cell with four iron sites. Such 
a match is the essential reason why the low energy physics 
becomes transparent after the gauge mapping. 

Third, the model unifies the iron-pnictides and iron- 
chalcogenides. When other parameters are fixed, reduc- 
ing t2.s or increasing tis can fiatten the dispersion along 
V — M direction of Eh± and cause the hole pocket com- 
pletely vanishes. Therefore, the model can describe both 
iron-pnictides and electron-overdoped iron-chalcogenides 
by varying t2s or tis- 

Finally, carefully examining the hopping parameters, 
we also find that the NNN hopping for each S'4 iso-spin 
essentially has a d-wave symmetry, namely \t2d\ > t2s- 
Since the hole pockets can be suppressed by reducing t2s 
value, this d-wave hopping symmetry is expected to be 
stronger in iron chalcogenides than in iron-pnictides. 

The two-orbital model with interactions: By projecting 
all interactions into these two effective orbital model, a 
general effective model that describes iron-based super- 
conductors obeying the S'4 symmetry can be written as 



Lq = 1.2 



Heff - Hq 

+C/'^n,4n,,2 + Jh^So -S.^s (11) 



where a = 1,2 labels the S4 iso-spin, U describes the 
effective Hubbard repulsion interaction within each com- 
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FIG. 5. A typical Fermi surfaces (a), band dispersions (b) re- 
sulted from Eq[8]with parameters t\s ~ 0.4, tid = —0.03, t23 = 
0.3, t2d = 0.6, t3s = 0.05, fad = -0.05 and ^J. = -0.3. (c) and 
(d) are corresponding results by adding tc = 0.02 in EqOwith 
the same parameter setting. 

ponent, U' describes the one between them and de- 
scribes the effective Hunds coupling. Since the two com- 
ponents couple weakly, we expect U dominates over U' 
and J'fi- In the first order approximation, the model 
becomes a single band-Hubbard model near half fill- 
ing. A similar t-J model can also be discussed within 
the same context as cuprates^^^. It is clear that the 
model naturally provides an explanation for the stable 
NNN AF exchange couplings J2 observed by neutron 
scattering^^^^^ and its dominating role in both mag- 
netism and superconductivity'^'^. 

Couplings between dxz,yz o,nd other d orbitals: Al- 
though the low energy physics is governed by dxz,yz or- 
bitals, it is interesting to ask their couplings to other d 
orbitals, in particular, the dxy orbitals. Under the S4 
transformation, the dxy is in the B representation while 
dxz,dyz are in the E representations of the S'4 group. 
Therefore, the coupling between them is only possible 
through other E representations, {sinkx, sinky) natu- 
rally forms two E representations. The leading coupling 
between dxz,yz and dxy is thus given by 

^xy — ix ^ ^ id^^(sinkxCk^a sinkydk^a} 

kcr 

+ty ^ idl^^isinkyCk.a + sinkxdk.a) + h.c, (12) 

ka 

where dka is the Fermionic operator for the dxy orbital. 
This form is uniquely determined by the time reversal 
symmetry and the ^4 symmetry. It is interesting to point 
out the coupling has been obtained where the authors 
tried to construct a three orbital modcl^ . This coupling 
can correctly account for the third hole pocket around 
(7r,7r) observed in the 122 structure^". The coupling to 



other d-orbitals can also be constructed similarly. In 
principle, the entire band structure can be constructed 
within the S'4 symmetry if necessary. 

Comparison between the S4 and the C^^ symmetry: It 
is interesting to compare the S'4 to the C^^ symmetry. If 
the Civ symmetry is forced, the C4 operation imposes, 

dk' + Qa \ (13) 
\dka J V -Ck' + Qo- J 

and the refiection operator in the 6*4^ group imposes 
additional requirements 

(t)-(r::)' <»' 

where k" — {ky, kx). It is easy to see that if we force the 
C4i, symmetry, the refiection invariance requires tis — 
0. However, without such a reflection invariance, this 
term is allowed, which is the case for the S'4 symmetry. 

The existence of tis suggests that (t„ symmetry must 
be broken in an effective model for d-orbitals. However, 
since symmetry appears to be present in Fe-As(Se) 
structure, it is natural to ask what mechanism can break 
(Jv While a detailed study of this symmetry breaking 
is in preparation^, we give a brief analysis. If we in- 
clude the (T„ together with the S4 group, the total group 
is enlarged to D2d- Therefore, in general, one can write 
down a D2d symmetry model by including all p orbitals of 
As(Se). When the model is reduced to an effective model 
based on the d-orbitals, because the p orbitals are fully 
occupied, the tis hopping between two NN dxzi dyz or- 
bitals can only be generated through the following virtual 
hopping processes: one electron first hops from the Px to 
the dxz, then, an electron in the Py at the same As(Se) 
site can hop to the Px, finally, an electron in the dyz or- 
bital hops to the Py. In such a process, the refiection 
symmetry is broken due to the existence of the hopping 
between the Px and py orbitals at the same As(Se) site 
when the two orbitals host total 3 electrons. 

The coupling between two S4 iso-spins and S'4 symme- 
try breaking: The couplings between the two iso-spins can 
either keep the S'4 symmetry or break it. Without break- 
ing the translational symmetry, the coupling between two 
orbitals can be written as 

Hc = Y.Uk)Go.{k) + J2M'')Gdk) (15) 

k.a k^a 

where Ga{k) and Ga{k) are operators constructed ac- 
cording to the S4 one dimensional representations as fol- 
lows, 

Gi (k) = 4adka + ct+Q^dk+Qa + h.c. (16) 

CT 

G2(fc) = Y4^dka - ct+Q^dk+Qa + h.c. (17) 

CT 

Gsik) = C+^dk+Qa + 4+Qa'^ka + h.C. (18) 
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FIG. 6. Fermi surfaces and band dispersions in the presence 
of the 5*4 symmetry breaking: (a, e) iti = 0.005 in Eq. 1241 
(b,f) tbt = 0.05 in Eq. [25l (c, g) Uo = 0.05 in Eq. [H (d,h) 
ihso = 0.05 in Eq. 1271 Other parameters are the same as in 
FigEl 





a 


h.c. 


(19) 


G-Ak) 


= ^^i4adka + 4+Qadk+Qa 
a 


- h.c.) 


(20) 


G-2{k) 


= '^ii4jka - 4+Qadk+Qa 
a 


- h.c.) 


(21) 


G-Ak) 


= H ^i4jk+Qa + 4+Qjka 
a 


- h.c.) 


(22) 


Giik) 


^J2^^4adk+Qa - 4+Qadka 


- h.c.) 


(23) 



We discuss a few examples that can cause the 5*4 sym- 
metry breaking, 

Hbi = ^ 2tbi{coskx + cosky){cl^dk+Qa + h.c.) (24) 

k 

Hbt AitbtsinkxSinky{c'l^dk+Qa - h.c.) (25) 

fc 

Hbo = ^ tboi4a^k+Qa - dta^k+Qcr) (26) 
fc 

Hbso ^^tbsoi4<7^ka - d^^dka)- (27) 
k 

tbi term breaks the S4 symmetry to Uft the degeneracy 
at r point, tf)t breaks the time reversal symmetry, tbo in- 
dicates a ferro-orbital ordering and tbso indicates a stag- 
gered orbital ordering. These terms can be generated 
either spontaneously or externally and their effects can 
be explicitly observed in the change of the band struc- 
ture and degeneracy lifting as shown in FiglHl where the 
changes of band structures and Fermi surfaces due to the 
symmetry breaking terms are plotted. It is fascinating to 
study how the interplay between the S4 symmetry and 
other broken symmetries in this system in future. 



IV. THE CLASSIFICATION OF THE 
SUPERCONDUCTING ORDERS ACCORDING 
TO THE 5*4 SYMMETRY 

The presence of the S4 symmetry brings a new symme- 
try classification of the superconducting phases. The S4 
point group has four one-dimensional representations, in- 
cluding A, B and 2E. In the A state, the 54 symmetry is 
maintained. In the B state, the state changes sign under 
the Si transformation. In the 2E state, the state obtains 
a ±7r/2 phase under the 6*4 transformation. Therefore, 
the 2E state breaks the C2 rotational symmetry as well 
as the time reversal symmetry. 

Since the S'4 transformation includes two parts, a 90° 
degree rotation and a reflection along c-axis, the S'4 sym- 
metry classification leads to a natural correlation between 
the rotation in a-b plane and c-axis reflection symmetries 
in a SC state. In the A-phase, rotation and c-axis reflec- 
tion can be both broken, while in the B-phase, one and 
only one of them can be broken. This correlation, in 
principle, may be observed by applying external symme- 
try breaking. For example, even in the A-phase and the 
rotation symmetry is not broken, we may force the c-axis 
phase-flip to obtain the phase change in the a-b plane. 

As shown in this paper, the iron-based superconduc- 
tors are rather unique with respect to the S'4 symmetry. 
It has two iso-spin components governed by the symme- 
try. This iso-spin degree of freedom and the interaction 
between them could lead to many novel phases. The fu- 
ture study can explore these possibilities. 

Here we specifically discuss the S4 symmetry aspects in 
the proposed Aig s-wave state, a most-likely phase if it is 
driven by the repulsive interaction or strong AF in iron- 
based superconductors'^?" as we have shown earlier. First, 
let's clarify the terminology issues. The Aig s-wave pair- 
ing symmetry is classified according to D^h point group. 
This classification is not right in the view of the true lat- 
tice symmetry. However, for each iso-spin components, 
we can still use it. Here we treat it as a state that the 
superconducting order A oc coskxcosk-^^. Since the Aig 
phase is equivalent to the d-wave in cuprates in a differ- 
ent gauge setting, the d-wave picture is more transparent 
regarding the sign change of the superconducting phase 
in the real space. As shown in FigdJ the sign of the SC 
order alternates between neighboring squares in the iron 
lattice. 

Based on the underlining electronic structure revealed 
here, with respect to the S4 symmetry, the Aig state can 
have two different phases, A phase and B phase. In the 
A phase, 

< Ck-fC^ki >=< dk\d_ki Aocoskxcosky (28) 
and in the B phase, 

< CfcfC-fc^ >= - < dk^d^ki >= Aocoshj:cosky. (29) 

Therefore, in the view of the d-wave picture, in both A 
and B phases, the superconducting phase for each com- 
ponent alternates between neighboring squares, which is 



r 




FIG. 7. (a) An illustration of a single Fe-As(Se) layer and 
the setup for a dc SQUIDS measurement to measure the sign 
change of the SC phase between top and down As(Se) layers; 
(b) The phase distribution in the A phase of the Aig S-wave 
state in the view of a d-wave picture ( red for one iso-spin 
component and blue for the other); (c) The phase distribution 
in the B phase of the Aig S-wave state. 



corresponding to the sign change between the top and 
bottom planes in the view of the S4 symmetry. However, 
in the A phase, since the 84 symmetry is not violated, 
the relative phase between the two components are equal 
to TT in space, while in the B phase, the relative phase is 
zero. A picture of the phase distribution of the two iso- 
spin components in the A and B phases are illustrated in 

Fig|7i;b,c). 

The sign change of the order parameter or the phase 
shift of TT between the top and bottom planes along c-axis 
can be detected by standard magnetic flux modulation 
of dc SQUIDS measurements^. If we consider a single 
Fe-As(Se) trilayer structure, which has been successfully 
grown by MBE recentlyi^ii^, we can design a standard 
dc SQUIDS as shown in FiglTlJa) following the similar 
experimental setup to determine the d-wave pairing in 
cuprates'S. For the B phase, there is no question that 
the design can repeat the previous results in cuprates. 
However, if the tunneling matrix elements to two com- 
ponents are not symmetric, even in the A phase, this de- 
sign can obtain the signal of the tt phase shift since the 
two components are weakly coupled and each of them 
has a TT phase shift. For the B phase, the phase shift may 
be preserved even in bulk materials^'i. However, for the 
A phase, it will be difficult to detect the phase shift in 
bulk materials. A cleverer design is needed. Measuring 
the phase shift between the upper and lower As planes 
will be a smoking-gun experiment to verify the model 
and determine iron-based superconductors and cuprates 
sharing identical superconducting mechanism. 
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FIG. 8. A sketch of the correlation between the hopping and 
pairing symmetries for both iron-based superconductors and 
cuprates. 



V. DISCUSSION AND SUMMARY 

We have shown that the Aig s-wave pairing in iron- 
based superconductors is a d-wave pairing in a view of 
a different gauge setting. This equivalence answers an 
essential question why a Aig s-wave pairing can be ro- 
bust regardless of the presence or absence of the hole 
pockets. With repulsive interactions, a sign changed or- 
der parameter in a superconducting state is usually in- 
evitable. This statement, however, is only true when the 
hopping parameters follow the same lattice symmetry. 
Gauge transformation can exchange the phases between 
superconducting order parameters and hopping parame- 
ters. In the case of cuprates, the d-wave order parameter 
can be transformed to a s-wave form by changing hopping 
parameters to obey d-wave symmetry. As we pointed out 
earlier, the NNN hopping in our model is close to a d- 
wave symmetry, rather than a s-wave symmetry. This is 
the essential reason why the superconducting order can 
have a s-wave form and be stable in iron-based super- 
conductors. A simple picture of this discussion is illus- 
trated in Fig[51 The vanishing of the hole pockets in 
electron-overdoped iron-chalcogenides indicates the hop- 
ping is even more d-wave like in these materials, a case 
supporting stronger s-wave pairing, which was indeed ob- 
served recentl y ^"'?!?. The presence of the dominant form 
coskxCosky is also straightly linked to the d-wave pairing 
form cosk'^ — cosk'y because of the stable AF J2 coupling, 
similar to cupates^. Moreover, since the different gauge 
setting does not alter physical measurements, a phase 
sensitive measurement should reveal a tt phase shift in 
the real space along c-axis for each components in the 
Aig s-wave state, just like the phase shift along a and b 
direction in the d-wave pairing state of cuprates. 

After obtaining the underlining electronic structure, 
we can ask how the physics in the cuprates and iron-based 
superconductors are related to each other. In TableU we 
list the close relations between two high Tc superconduc- 
tors. From the table, it is clear that by determining these 
physical properties of iron-based superconductors listed 
in the table can help to determine the high supercon- 
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properties 


Iron SCs 


cuprates 


pairing symmetry 


s-wave 


d-wave 


underlining hopping symmetry 


d-wave 


s-wave 


dominant pairing form 


coskxCosky 


coskx — cosky 


AF coupling 


NNN Ja 


NN Ji 


sign change in real space 


c-axis 


a-b plane 


filling density 


half-filling 


half-filling 



TABLE I. A list of the close connections between iron-based 
superconductors (iron SCs) and cuprates. 



ducting mechanism. 



The model completely changes the view of the origin 
of the generation of sign-changed pairing symmetry in 
iron-pnictides, which were argued in many theories that 
the origin is the scattering between electron pockets at 
M and hole pockets at F due to repulsive interactions^iS,. 
With the new underlining electronic structure revealed, 
the analysis of the sign-change should be examined after 
taking the gauge transformation so that the underlining 
hopping parameters become symmetric. In this case, the 
sign change is driven by scatterings between all pockets, 
including both hole and electron pockets, located at two 
d-wave anti-nodal X' points. Therefore, the scattering 
between electron and electron pockets is also important. 

While the model appears to be rotational invariant due 
to the S4 symmetry, the dynamics of each iso-spin compo- 
nent is intrinsically nematic. A small S'4 symmetry break- 



ing can easily lead to an overall electronic nematic state. 
The electronic nematic state has been observed by many 
experimental techniques^ and studied by different theo- 
retical models^"—. The underlining electronic structure 
in the model can provide a straightforward microscopic 
understanding between the interplay of all different de- 
gree of freedoms based on the S'4 symmetry breaking. 

In summary, we have shown the underlining electronic 
structure, which is responsible for superconductivity at 
low energy in iron-based superconductors, is essentially 
two weakly-coupled electronic structures governed by the 
^4 symmetry. We demonstrate the s — wave pairing in 
iron-based superconductors is equivalent to the d-wave 
in cuprates. A similar conclusion has also been reached 
in the study of 2-layer Hubbard model?'*. The ^4 sym- 
metry reveals possible new superconducting states and 
suggests the phase shift in the SC state in real space 
is along c-axis. These results strongly support the mi- 
croscopic superconducting mechanism for cuprates and 
iron-based superconductors are identical, including both 
iron-pnictides and iron-chalcogenides. Our model estab- 
lishes a new foundation for understanding and exploring 
properties of iron based superconductors, a unique, ele- 
gant and beautiful class of superconductors. 
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